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ABSTRACT 

Deep XMM-Newton and Chandra observations of C1J1226. 9+3332 at z = 0.89 have enabled the most 
detailed X-ray mass analysis of any such high-redshift galaxy cluster. The XMM-Newton temperature 
profile of the system shows no sign of central cooling, with a hot core and a radially declining profile. A 
temperature map shows asymmetry with a hot region that appears to be associated with a subclump 
of galaxies at the cluster redshift, but is not visible in the X-ray surface brightness. This is likely 
to be result of a merger event in the cluster, but does not appear to significantly affect the overall 
temperature profile. The XMM-Newton temperature profile, and combined Chandra and XMM- 
Newton emissivity profile allowed precise measurements of the global properties of C1J1226. 9+3332; 
wc find kT = 10.4 ± 0.6 keV, Z = 0.16 ± O.O5Z , and M = 5.2+*;° x 10 14 M o . We obtain profiles of 
the metallicity, entropy, cooling time and gas fraction, and find a high concentration parameter for 
the total density profile of the system. The global properties are compared with the local L — T and 
M — T relations, and we are able to make the first observational test of the predicted evolution of the 
Yx — M500 relation. We find that departures from these scaling relations are most likely caused by an 
underestimate of the total mass by ~ 30% in the X-ray hydrostatic mass analysis due to the apparent 
recent or ongoing merger activity. 

Subject headings: cosmology: observations - galaxies: clusters: individual: (C1J1226. 9+3332) - galax- 
ies: clusters: general - galaxies: high-redshift galaxies: clusters - intergalactic 
medium - X-rays: galaxies 



1. INTRODUCTION 

Clusters of galaxies are an important tool for test- 
ing and improving cosmological models. High-redshift 
clusters, in particular, provide the strongest constraints. 
Measurements of the total and baryonic mass functions, 
and baryon fractions of distant clusters have been used 
to measure cosmolo gical parameters (e .q.lVikhlini n et"aT1 
120031 IHenrvl l200l lAllen et alJ l200l . Even with the 
precise measurements of cosmological parameters pro- 
vided by cosmic microwave background and supernov a 
observations fe.o. lSpergel et al1l2003t iTonrv et aLl feOOS 1 ). 
clusters play an important role by providing indepen- 
dent constraints with different parameter degeneracies 
tl Allen et al.ll2004|) . Large samples of distant clusters de- 
tected in planned X-ray and Sunyaev-Zel'dovich effect 
(SZE) surveys will be able to place strong constraints 
on the dark energy density and equation of state, but 
the most powe rful constraints require ma ss estimates for 
many clusters l)Maiumdar fc Mohrll2003ft . 

The most reliable X-ray mass estimates require ra- 
dial temperature and density profiles of the intra-cluster 
gas to allow the solution of the equation of hydrostatic 
equilibrium. However, the low X-ray fluxes of typical 
high-redshift clusters make the observing time needed to 
meet these requirements prohibitive. Sufficiently deep 
observations for clusters at z > 0.5 are very rare (e.g. 
[ Jeltema et all 120011 lArnaud et alJl2002t iDonahue et all 
120031 iMauehan et all l2004f ~ for most clusters at high 
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rcdshifts, masses must be estimated from observed global 
properties such as X-ray luminosity and temperature. 
The scaling relations between cluster properties at high 
rcdshifts are then of key importance, and their uncer- 
tainty is a dominant contributor to the error budget of 
derived cosmological parameters ijHenrvl 12004) . Further- 
more, studies of the evolution of the cluster scaling rela- 
tions can be used to probe cluster for mation and the his- 
tory of non-gravitatio nal heating fe.o. lEttori et al.ll2003t 
iMaughan etaI1l2006[) . 

At z = 0.89, galaxy cluster C1J1226.9+3332 is the 
hottest, most luminous known at z > 0.6. 4 The only 
comparable system is MS1054.5— 0321 at z — 0.83, an ac- 
tively merging system which is s omewhat cooler and less 
luminous than C1J1226.9+3332 (|Jee et al.ll2005|) . Unlike 
that very disturbed system, C1J1226. 9+3332 has a re- 
laxed morphology, making it a promising candidate for 
a reliable hydrostatic mass estimate. C1J1226. 9+3332 
was detec ted in the Wide Angle ROSAT Pointed Survey 
fWARPS lScharf et aljfl997t lEbeling : et alJl2001|) and an 
analysis of a 15 ks AOl XMM-Newton observation con- 
firmed its temperature at 11.5 keV with a bolometric lu- 
minosity of 5 x 10 45 erg s _1 and a relaxed morphology 
l)Maughan et alJ2004[) . An early, short Ch andra observa- 
tion o f C1J1226. 9+3332 was presented bv lCagnoni et all 
1)200 1|) . who also pointed out that the VLA FIRST cat- 
alogue shows a faint point-like radio source at the clus- 
ter centre, consistent with a low-luminosity radio-loud 
AGN C1J1226 9+3332 has also been detected in the 
SZE l)Jov et al.ll2001[) and the properties of the mem- 

4 Comparisons here are made with all clusters listed in the 
BAX database I http : //bax . ast . obs-mip. f r I, a compilation of 
published X-ray properties of galaxy clusters 
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ber galaxies hav e been investigated ijEllis &: Jonesll200"l 
Ellis et alJl200fil) . 

In this paper we present new XMM-Newton and Chan- 
dra observations of C1J1226. 9+3332. In the following 
sections we describe the data reduction and background 
subtraction techniques, and present the results of the 
imaging and spectral analysis. The radial profiles of the 
clusters properties such as mass, temperature and en- 
tropy are then investigated and the global properties of 
the system are compared with those of other clusters. 
In the appendix we present a novel method for defin- 
ing contour levels in smoothed X-ray images. Numerical 
subscripts are used on cluster properties to indicate the 
overdensity with respect to the critical density at the 
radius at which (or within which) the property was mea- 
sured. The critical density at the cluster's redshift is used 
in this definition. We assume a flat ACDM cosmology 
with H = 70 km s" 1 Mpc -1 and Q M — 0.3 throughout, 
and uncertainties are quoted at the 68% level. 

2. DATA REDUCTION 

The XMM-Newton data (ObsID 0200340I0I) were re- 
duced and analysed with the Science Analysis Software 
(SAS) version 6.5, and the most recent calibration prod- 
ucts available as of 15 November 3005. A lightcurve of 
the field of view was produced in the 10 — 15 keV band 
for each of the three EPIC cameras. The lightcurves 
were filtered to reject periods of high background, and 
approximately 16 ks of data were rejected. The cleaned 
lightcurves exhibited stable mean count rates and the 
clean exposure times were 68 ks and 75 ks for the PN and 
MOS de tectors respectively. A l ocal background region, 
and the iRead fc PonmarJ l)2003fl blank-sky background 
datasets were used for background estimates where ap- 
propriate, as discussed below. 

The Chandra data consisted of two ACIS-I observa- 
tions (ObsID 3180 and 5014) which were were reduced 
and analysed consistently. The level 1 events files were 
reprocessed with the latest calibration (as of November 
2005), applying the charge-transfer inefficiency and time- 
dependent gain corrections. Background flare filtering 
was performed with lc_clean in the 0.3 — 12 keV band 
excluding the target CCD, and all sources detected by 
the pipeline processing. The good time remaining for 
each observation was 25 ks. The period D blank-sky 
background files l)Markevitchll200 3) were appropriate for 
both observations, and these were reprojected onto the 
sky coordinates of each observation. Both observations 
were taken in VFAINT telemetry mode, so the additional 
VFAINT cleaning procedure was applied to the source 
and background datasets 5 . 

3. BACKGROUND CONSIDERATIONS 

In periods free of background flaring, the background 
is dominated by X-rays at low energies ( < 2 keV) and 
particles at high energies ( > 2 keV). The X-ray com- 
ponent includes a significant contribution from Galac- 
tic emission which varies with position on the sky. The 
blank-sky background datasets are compiled from mul- 
tiple pointings with different Galactic X-ray fluxes, so 
may not agree well with the soft X-ray background in the 
source dataset. However, the blank-sky datasets enable 

5 http:/ /asc. harvard.edu/ciao/threads/aciscleanvf 
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Fig. 1. — Source and background radial surface brightness profiles 
of C1J1226. 9+3332 measured with XMM-Newton and Chandra. 



background estimates to be made at the same detector 
position as the source, eliminating instrumental effects, 
and furthermore, many extended sources fill the field of 
view of Chandra or XMM-Newton prohibiting the mea- 
surement of the background in the source dataset. Due 
to its high redshift, C1J1226. 9+3332 occupies a small 
area on the detectors, leaving large source-free regions in 
which the background emission can be measured. This 
local background was compared with the blank-sky files 
to determine the optimal choice of background for imag- 
ing and spectral analyses. 

Images were created in the 0.7 — 2 keV band for 
each detector, along with exposure maps to correct for 
chip gaps, bad pixels, and telescope vignetting. The 
three XMM-Newton PN and MOS images and the two 
Chandra ACIS-I images were combined to give sepa- 
rate XMM-Newton and Chandra images, and these were 
exposure corrected. Background images and exposure 
maps were created in the same way from the blank- 
sky datasets for each detector. Sources in these images 
w ere then detected with t he wavelet detection algorithm 
of iVikhlinin et al. (199!§). The detected point and ex- 
tended sources were then excluded from the analysis of 
C1J1226.9+3332. 

Due to the possible difference in the soft Galactic X- 
ray flux between the source and background datasets, the 
blank-sky images were normalised to match the flux in 
the source-free regions in the target images. For both the 
XMM-Newton and Chandra images, radial profiles of the 
source and background images and their exposure maps 
were created and used to determine the detection radius 
(rd) of the cluster emission. Emission within rd was then 
excluded and the background images were normalised to 
the remaining emission. The detection radius was the 
recomputed and the process was repeated until rd con- 
verged. This procedure provides a reliable measurement 
of both the cluster extent and the normalisation of the 
blank-sky images in the imaging band. 

Fig. ^ shows the source and background profiles for 
the XMM-Newton and Chandra images. The good agree- 
ment of the source and background profiles at r > 300" 
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indicates that the source exclusion and background nor- 
malisation were successful. Note that the increasing 
background with radius (most noticeable in the XMM- 
Newton data) is due to the fact that the particle back- 
ground component is not vignetted, so the exposure cor- 
rection incorrectly boosts this background component. 
In this imaging analysis, the effect is the same in the 
source and background images and cancels out. How- 
ever, this effect is a cause for concern if an off-axis region 
is used to extract a local background spectrum. These 
normalised blank-sky background images were thus used 
for all of the imaging analysis for all datasets, while the 
background subtraction for the spectral analysis is dis- 
cussed in more detail below. The consistency of the 
XMM-Newton and Chandra surface brightness profiles 
of C1J1226. 9+3332 is demonstrated in section iJ6.ll 

A prime motivation in defining rd was to fully exclude 
any significant cluster emission to enable to measurement 
of the background in the target observations. To this end 
rd was conservatively defined as the radius beyond which 
no further radial bins had a detection of > 0.5<r above 
the background. The detection radius for XMM-Newton 
was 210" and for Chandra was 126". For each detector, 
spectra were extracted from the target datasets after ex- 
cluding all sources, and the cluster emission out to rd, 
to give a local background spectrum. Blank sky back- 
ground spectra were extracted for each detector with the 
same spatial regions excluded, and the local and blank 
sky background spectra are plotted in Fig. H 

The Chandra blank-sky background spectra were nor- 
malised to match the 9.5 — 12 keV count rates in the 
target datasets. At these high energies, the background 
is dominated by particle events. The left panel of Fig. 
12 shows the background spectra for the Chandra Ob- 
sID 3180 observation (ObsID 5014 is consistent). The 
good agreement between the spectra above 1.5 keV indi- 
cates that relative contributions of the particle and X-ray 
background components is the same in the target and 
blank-sky datasets. The excess in the local background 
at the soft end is due to Galactic X-ray emission. These 
soft residuals are accounted for in our spectral analy- 
sis of the Chandra data by including a model for this 
soft thermal component (see SJSJ. We conclude that the 
blank-sky background files provide a reliable background 
measurement for the Chandra datasets. 

The XMM-Newton blank-sky background spectra were 
normalised by the ratio of the count rates of events de- 
tected outside the field of view of the XMM-Newton 
telescopes. These events are due solely to the particle 
background. It is clear in the centre and right panels 
of Fig. [21 that the flux in the particle-induced fluores- 
cent lines agrees well in the target and blank-sky spec- 
tra, but the flux elsewhere is significantly lower in the 
local background spectra. This suggests that the rela- 
tive contribution of the particle background component 
is significantly greater in the target dataset. A similar ef- 
fect has been found in other XMM-Newton obs ervations 
ijMaughan et alJl2006t iKhosroshahi et al.ll2006j) . and we 
emphasise the importance of comparing the blank-sky 
and local background spectra whenever possible. A pos- 
sible resolution to this discrepancy would be to normalise 
the spectra by the count rate in the 2 — 7 keV band, but 
this would lead to an underestimate of the flux in the 
fluorescent lines. Furthermore, the slopes of the of the 



spectra are different due to the differing particle back- 
ground contributions. We thus conclude that the blank- 
sky datasets cannot be used for the spectral analysis of 
this XMM-Newton observation. Instead, we define a lo- 
cal background annulus region between 250 < r < 500". 

This local background region lies further off-axis than 
the cluster emission, so the vignetting of the background 
photons must be considered. Vignetting effects were cor- 
rected by using the SAS task evigweight. This assigns a 
weight to every event based on the ratio of the effective 
area at its detected position and energy, to the effective 
area on-axis at the same energy. This procedure can not 
distinguish between vignetted X-rays, and non- vignetted 
particle events, leading to the incorrect boosting of the 
particle background. Fig. ^ demonstrates that this is 
already a significant effect by r ~ 400". 

For a given source region, this effect can be removed 
with the aid of the XMM-Newton closed filter back- 
ground datasets which consist of purely instrumental and 
particle events. A "boosted-particle" spectrum was cre- 
ated for each camera by running evigweight on the closed 
filter background and extracting a spectrum from the lo- 
cal background region. This replicates the treatment of 
the particle events in the source dataset. These boosted- 
particle spectra were then normalised to the 10 — 14 keV 
count rate in the local background spectra and sub- 
tracted from them, creating a vignetting-corrected "X- 
ray only" background spectrum. A similar boosted- 
particle spectrum was then extracted from the source 
region in the closed filter data, normalised by the same 
factor and added onto the "X-ray only" background spec- 
trum (after scaling by area). The resultant background 
spectrum then consisted of an X-ray background that 
was correctly weighted for the vignetting at the back- 
ground region, and a particle component that was in- 
correctly weighted in exactly the same way as the par- 
ticles in the source spectrum. The net effect is that the 
particle and X-ray backgrounds are both subtracted cor- 
rectly. This has a significant effect on the temperatures 
measured at large radii where the background dominates, 
and this method was followed for all XMM-Newton spec- 
tral analysis. 

4. CLUSTER MORPHOLOGY 

The exposure corrected Chandra and XMM-Newton 
images of C1J1226. 9+3332 were smoothed to the 3cr level 
using the adaptive smoothing algorithm of Eb eling et alJ 
(200$). Contours were then derived from the smoothed 
images at flux levels defined so that the emission in the 
region between two adjacent contour levels was detected 
at > 3cr above the emission in the next lowest band. This 
contouring algorithm is described in detail in Appendix 
Eland allows the significance of features in contour plots 
to be determined easily by eye. 

Fig. Olshows a wide field view of C1J1226.9+3332 with 
XMM-Newton contours overlaid on a Keck-II I band 
image. The image is dominated by the bright, circu- 
lar C1J1226. 9+3332 in the western half, but the fore- 
ground cluster C1J1227. 3+3333 is visible in the eastern 
half of the image. C1J1227. 3+3333 appears unrelaxed 
and consists of several sub-clumps of galaxies at a simi- 
lar spectroscopic redshift (z = 0.766). All of the sources 
seen in the contours are consistent with point sources 
except for the northern peak of C1J1227. 3+3333 and 
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Fig. 2. — Background spectra taken from the target (black points) and blank-sky (grey points) datasets for the Chandra ACIS-I ObsID 
3180 (left), XMM-Newton PN (centre) and XMM-Newton MOS (right) observations of C1J1226. 9+3332. 
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Fig. 3. — Contours of the emission detected by XMM-Newton 
are overlaid on a Keck-II I band image. The foreground cluster 
C1J1227. 3+3333 (z = 0.766) is visible in the eastern part of the 
image. 




777 kpc 



08.0 06.0 04.0 02.0 12:27:00.0 50.0 56.0 54.0 52.0 26:50.0 45.0 
Right ascension 

Fig. 4. — Contours of the emission detected by Chandra are 
overlaid on a Keck-II I band image. 



C1J1226. 9+3332. Fig. H shows the Chandra contours 
of C1J1226. 9+3332 on the same optical image. These 
contours demonstrate the regular, relaxed appearance of 
the emission in C1J1226.9+3332. 

5. SPECTRAL ANALYSIS 
5.1. Methodology 



Before discussing the different spectral analyses per- 
formed, we give an overview of the spectral extraction 
and fitting methods used for the XMM-Newton and 
Chandra data. In all cases, responses were generated 
for the spectral extraction regions, weighted by the dis- 
tribution of counts in the region. Cluster spectra were 
modelled in XSPE C with a single temperature, APEC 
(Smit het alJl20 01') plasma model, and with the excep- 
tion of the MOS data as noted below, the absorbing col- 
umn was fixe d at the Galactic value inf erred from 21 cm 
observations l)Dickev fc LockmarJll990l) . During fitting, 
the APEC normalisation, metal abundance and temper- 
ature were the only free parameters. When data from 
multiple detectors were fit simultaneously, the tempera- 
ture and metal abundances of the models were tied, but 
the normalisations were not. 

For the XMM-Newton data, only the best-calibrated 
events (FLAG=0) were selected, and the modified local 
background spectra discussed in fJ3I were used for each 
detector. The MOS and PN spectra were fit in the 0.4 - 
7 keV energy bands, avoiding the bright fluorescent lines 
above 7 keV in the PN and minimising any remaining 
problems with the particle background subtraction. 

For the Chandra data, the background spectra were 
taken from the blank sky datasets. A spectrum of the 
soft residuals due to Galactic X-ray emission differences 
between the target and blank-sky datasets was produced 
for each dataset by subtracting a blank-sky background 
spectrum from a local background spectrum extracted 
in the same region. This was fit with a 0.18 keV, un- 
absor bed APEC model of zero redshift (Vikhlini n et alJ 
2005). This model was included as an additional back- 
ground component in all spectral fits, with the normal- 
isation scaled by the difference in extraction area. The 
cluster spectrum was fit by an absorbed APEC model 
in the 0.6 — 9 keV and the best fitting temperature was 
found. The model was refit with the soft Galactic compo- 
nent normalisation set to ±lcr and again with the overall 
blank-sky background renormalised by ±2%, and the re- 
sulting systematic temperature uncertainties were added 
in quadrature to the statistical uncertainties. 

5.2. Cross- calibration 

In order to test the spectral calibration of the different 
detectors before combining their data, a single spectrum 
was extracted from r < 50" for each detector. Select- 
ing only the central cluster region in this way maximised 
the signal to noise ratio of the spectra, reducing the un- 
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Fig. 5. — Spectra extracted from the central 50" of 
C1J1226. 9+3332 plotted with best fitting models with Z = O.3Z . 
Top: XMM-Newton spectra; the temperature of the 3 models were 
tied together, and the obsorbing column was fixed at the Galac- 
tic value for the PN data and at 4.1 X 10 20 cm" 2 for the MOS 
data (lower lines). Bottom: Chandra spectra with soft residuals 
modelled simultaneously. The temperatures of the Chandra cluster 
spectra were tied together. 



certainties due to the background subtraction. No point 
sources were detected within this region by either Chan- 
dra or XMM-Newton. To further simplify the compar- 
ison, the metal abundances were fixed at 0.3Z©. The 
spectra were fit separately and simultaneously in differ- 
ent combinations, and example spectra are plotted in 

Fig. El 

The results of the spectral fitting are summarised in 
Fig. HO The solid points show that the temperatures 
measured with EPIC MOS were significantly higher than 
the PN measurements; the 90% confidence intervals do 
not overlap. The absorbing column was then allowed 
to vary, and the agreement in kT was improved dramati- 
cally. The best fitting PN absorbing column did not vary 
significantly from the Galactic value (1.38 x 10 20 cm" 2 ) 
and the PN temperature was unchanged. In the case of 
the combined MOS spectra, however, the best-fitting ab- 
sorption was significantly higher than the Galactic value 
at 4.1±1.0xl0 20 cm" 2 (90% errors) and the temperature 
was significantly reduced. The absorbing columns found 
when fitting the MOS data independently of each other 
were also consistent with this value. All of these spectral 
fits were statistically acceptable, whether the absorbing 
column was fixed or varying. 

The hollow points in Fig. [B] show the tempera- 
tures measured when the absorbing column was fixed at 
4.1 x 10 20 cm" 2 for the MOS data; all EPIC measure- 
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Fig. 6. — Comparisons of temperatures measured for the central 
50" of C1J1226. 9+3332 with XMM-Newton and Chandra detec- 
tors. The labels indicate the detector combination and confidence 
level for each measurement. The XMM-Newton spectra were fit in 
the 0.4 — 7 keV band. Solid points indicate temperatures measured 
with the absorbing column fixed at the Galactic value; for the hol- 
low points the absorbing column for the MOS detectors was fixed 
at 4.1 x 10 20 cm" 2 . 



ments agree well. The statistical precision of the ACIS 
measurements, meanwhile, cannot exclude any of the 
XMM-Newton measurements. The fact that the MOS 
fits require an absorbing column significantly higher than 
Galactic suggests that there is an outstanding calibration 
problem with those detectors, and that the PN tempera- 
ture should be more reliable. In all further spectral fits, 
the absorbing column was frozen at the Galactic value 
for the PN and at 4.1 x 10 20 cm" 2 for the MOS data. 
The variation in the best-fitting temperature caused by 
setting the MOS absorbing column to its 90% limits was 
taken as an indication of the additional systematic uncer- 
tainty due to this problem. This uncertainty was added 
in quadrature to the statistical errors on all subsequent 
temperature measurements, resulting in an increase of 
w 50%. Finally, we note that consistent temperatures 
and uncertainties were found using the combined EPIC 
data if energies below 1.2 keV were ignored in the MOS 
spectra. 

5.3. Spectral Profiles 

The XMM-Newton observation of C1J1226. 9+3332 is 
sufficiently deep to provide radial profiles of its spectral 
properties. Spectra were extracted from annuli centred 
on the X-ray centroid (with point sources excluded) , and 
were fit as before. The resulting projected temperature 
profile is plotted in Fig. [7| The data from the three EPIC 
camera were fit simultaneously for these results, with the 
temperatures tied, the PN absorbing column frozen at 
the Galactic value, and the MOS absobing column at 
4.1 x 10 20 cm" 2 

To model the tempera ture profile of C1J1226 . 9+3332 
we follow the method of IVikhlinin et alJ (|2006l V06) of 
projecting a model of the three-dimensional (3D) tem- 
perature distribution along the line of sight, and fitting 
that to the observed projected tem perature profile. In 
doing this we use the algorithm of Vikn lmhl (f2006^ to 
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Fig. 7.— The projected temperature profile of C1J1226.9+3332 
measured with XMM-Newton is plotted along with the best fitting 
3D model and its projection. 



Radius (kpc) 

Fig. 8.— Projected metal abundance profile of C1J1226.9+3332. 
The y errorbars show the 68% confidence intervals. 



predict the observed temperature for the combination 
of temperatures and emissivities of gas along the line 
of sight in each bin. The emissivities were determined 
from our analysis of the gas density profile (JOJ, and 
the metal abundance of the gas was set at O.3Z at all 
radii (varying the metal abundance had negligible effect 
on this analysis, which is reasonable, as the emission is 
continuum dominated at these high temperatures). The 
temperature profile model of V06 was used to model the 
3D temperature distribution: 

kT(r) _ kx (r/r cool r^+kT mm /kT Q (r/r t )~° 

This model was simplified slightly to reduce the number 
of free parameters by fixing a = and b = c/0.45, in 
line with the average profile found by V06 in a sample 
of relaxed nearby clusters. The best fitting 3D model 
had kT = 22.0 keV, r cool = 47.4 kpc, a cool = 0.0, 
kT min = 2.7 keV, r t = 329 kpc, and c = 0.72, and the 
model and its projection are plotted in Fig. This 
model has as many degrees of freedom as datapoints in 
the profile, and has degeneracies between parameters, so 
the uncertainties on all dependent quantities were com- 
puted by a Monte-Carlo method. A simpler model 



kT(r)=kT (l + r/r t y 



(2) 



was also fit to the data, giving kTo = 13.3 keV, r t = 
1.0 Mpc and a = 1.1. While the V06 model was used for 
all results quoted subsequently, this simpler model gave 
consistent results and uncertainties throughout. 

This method of modelling the temperature profile 
does not take the effect of the XMM-Newton PSF 
into account. As the surface brightness profile of 
C1J1226. 9+3332 is not sharply peaked and the tempera- 
ture profile bins are all > 15" (the approximate FWHM 
of the PSF), the mixing effect of the PSF on the tempera- 
ture profile is not expected to be large. This was verified 
by fitting the spectra from all radial bins simultaneously 
with absorbed APEC models, using the xmmpsf model 
in XSPEC to model the redistribution of photons by the 



XMM-Newton PSF. The resulting PSF-corrected tem- 
perature profile was consistent with the original profile, 
with the noise fluctuations amplified and large uncertain- 
ties. The original uncorrected temperature profile was 
used for all further analysis. 

The projected metal abundance profile of 
C1J1226.9+3332 is plotted in Fig i] The statisti- 
cal uncertainties on the metal abundance measurements 
are large due to the high ionisation state of the plasma. 
The data suggest a radially decreasing abundance profile 
but are also consistent with a flat abundance profile. 

5.4. Temperature map 

The two dimensional projected temperature distribu- 
tion in C1J1226. 9+3332 was investigated by constructing 
a temperature map fo llowing; the methods described in 
iMaughan et~aT1 (|2006t see also IP 'Sullivan et alJ <|2005[n . 
In brief, spectra were extracted from many overlapping 
circular regions in the cluster, the radii of which were 
chosen to contain at least 1000 net photons. The best 
fitting temperature for the spectrum from each region 
was assigned to the pixel in the temperature map cor- 
responding to the centre of that region. The result is a 
smooth image of the projected temperatures, in which 
nearby pixels are not independent. The spectra were fit 
as before, but the metal abundance was fixed at 0.3Z Q . 

The resulting projected temperature map is shown in 
Fig. |5J The temperature structure shows asymmetry 
that is not apparent in the surface brightness distribu- 
tion. There is significantly hotter emission to the west 
and south west of the cluster core than the surrounding 
gas at the same radius. Fig. ^] shows contours of the 
X-ray temperature overlaid on the optical image. The re- 
gion of hotter emission appears associated with an over- 
density of galaxies (confirmed cluster members) to the 
south west of the core, suggesting that a merger event 
may be responsible for the temperature asymmetry. 

A new temperature profile was produced, with the sec- 
tor containing these hotter regions excluded, and the 
temperatures obtained were not significantly different 
from those measured with the sector included. This in- 
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dicates that this hot region does not significantly bias 
our temperature profile measurement. However, if it is 
related to merger activity in the cluster, this will under- 
mine the assumption of hydrostatic equilibrium used in 
the following mass analysis of C1J1226. 9+3332. 

6. MASS ANALYSIS 
6.1. Gas density profile 

The X-ray emissivity of the intra-cluster medium de- 
pends strongly on the gas density and only weakly on its 
temperature. This means that the observed projected 
emissivity profile can be used to accurately measure the 
gas density profile. For each annular bin in the surface 
brightness profile, the observed net count rate was cor- 
rected for area lost to chip gaps, bad pixels and excluded 
regions, and converted to an integrated emission mea- 
sure. For each bin, the conversi on factor was calculated 
assuming an absorbed MeKaL ijKaastra fc Mewdll993]) 
plasma model folded through an ARF generated for that 
region, and an on-axis RMF. The absorption of the spec- 
tral model was set at the galactic value and the metal 
abundance was set at 0.3. The temperature of the spec- 
tral model in each bin was given by the best-fitting 3D 
temperature model. We note that the dependence of the 
conversion from count rate to emission measure on the 
assumed temperature is weak; this analysis was repeated 
assuming the cluster was isothermal for the purpose of 
deriving the emission measures, and the results were not 
significantly changed. 

Projected emission measure profiles were derived from 
the ACIS and PN data in this way, and are plotted in Fig. 
1111 The radial bins were set to a minimum width of 15" 
for the PN profile, to minimise the effect of the XMM- 
Newton PSF, and the profiles agree very well. The MOS 
data were omitted here for simplicity, as the PN data 
alone were sufficient to extend the profile beyond the 
range of the temperature profile, which is the limiting 
factor in the mass analysis. 

To model the emissivity profile, we used the modified 
/3-profile used by V06 to describe the emission measure 
profile of the gas in nearby relaxed clusters; 

We = nl (1 + ( ;/ / r ;i ) 7,- a/2 (1 + rVr])-^. (3) 

This model can describe a steep core component and any 
steepening of the data beyond the standard /3-profile at 
large radii. A additive core component has been removed 
from the V06 model here as it was not required by the 
data. This emission measure model was then projected 
along the line of sight and fit to the observed projected 
emission measure profile. The model was fit to the Chan- 
dra and XMM-Newton PN data simultaneously, with the 
PN data within 500 kpc excluded to minimise PSF ef- 
fects. The best fitting model is plotted in Fig. The 
gas density is then related to the emission measure by 
p g = 1.252m p (n p n e ) 1 / 2 , assu ming a cosmic plasma wit h 
helium abundances given bv lAnders fc Grevessd l)1989|) . 
For comparison with other work, the profile was also fit 
with a standard /3-model, and the best-fitting parameters 
were r c = 105 ± 5 kpc and f3 = 0.64 ± 0.01. 

6.2. Mass profiles 



The 3D models of the gas density and tempera- 
ture distributions enable the total gravitating mass of 
C1J1226. 9+3332 to be derived under the assumptions of 
spherical symmetry and hydrostatic equilibrium. Fig. 1121 
shows the total and gas mass profiles of C1J1226. 9+3332 
with the shaded regions indicating the la uncertainties. 
The error analysis was performed by producing 1000 ran- 
domisations of the observed temperature and emission 
measure profiles based on their measurement uncertain- 
ties. The full analysis was repeated for each randomisa- 
tion, and the ±34 percentiles of the distribution about 
the value determined from the non-randomised data gave 
the ±lcr uncertainties on each parameter. 

The enclosed gas mass fraction (/ gas ) profile of 
C1J1226. 9+3332 was also computed, and the resulting 
profile is plotted in Fig. The profile increases out 

to -R500, which is consistent with the profiles measured 
b y V06 for local clus ters, but contrasts with the findings 
of I Allen et al. (200J) that / gas profiles tend to a univer- 
sal value at around -R2500 (although with some variation 
between clusters). The latter work used Chandra obser- 
vations of sample of clusters including C1J1226. 9+3332; 
the value of / gas we find for C1J1 226. 9+3332 at -B2R 00 is 
consistent with that measured by I Allen et al.l l)2004t) . 

7. ENTROPY AND COOLING TIME PROFILES 

The entropy (K) of the intra-cluster gas can give 
usefu l insight into its ener getics (e.g. iPonman et alJ 
119991 iVoitfc Donah ue 200|). The entropy profile of 
C1J1226. 9+3332 was derived from the 3D temperature 
and gas density models according to the definition K = 

kTn e 2 ^ 3 and is plotted in Fig. In order to model 
this entropy profile, we started with the baseline en- 
tropy profile due to gravitati onal effects alone , derived 
from numerical simulations bv lVoit et alJ l)2005[k K(r) — 
1.32^200 {f / R200) 1 ' 1 • K-ion was estima ted following the 
method suggested bv lVoit et al.l l)2005|) . with 

if 2 oo = 362fc!r 2 oo£(z)- 4 / 3 , (4) 

where E(z) gives the redshift-variation of the Hub- 
ble parameter and fcX^oo was estimated from the mass 
within R200 via fcT 20 o = GM 2 oo^m p /2R2oo. For 
C1J1226.9+3332 kT 2Q0 was found to be 6.6 keV by ex- 
trapolating the mass profile model to R 2 00i aud hence 
K 2 qq = f226 keV cm 2 . The baseline entropy profile is 
plotted in Fig. Upland is clearly not a good description 
of the data. The data suggest a large entropy excess 
above the gravitational baseline in the central regions of 
C1J1226. 9+3332. To model this, a constant entropy level 
(Kq) was added to the baseline model to give 

K(r) = K + 1.32^2oo(r/i?2oo) 1 ' 1 . (5) 

This model was fit to the entropy values mea- 
sured at the midpoints of the 6 temperature bins of 
C1J1226.9+3332. The best fitting model had K = 
132 ± 24 keV cm 2 and is plotted in Fig. d The 
model is a good description of the data (x 2 jv — 1-7/5) 
and a s ignificant impro v ement on the baseline model 
alone. iDonahue et alJ (2006) used a similar model 
to fit entropy profiles in nearby cooling flow clusters, 
and the central entropy level Kq is significantly higher 
in C1J122 6.9+3332 than the values of ~ 10 keV cm 2 
found bv IDonahue et afl l|2006j) . This is consistent with 
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Fig. 9. — Left: Projected temperature map of C1J1226. 9+3332 with Chandra contours of X-ray surface brightness overlaid. Right: The 
upper panels show the lc upper (left) and lower (right) limits on the temperature values. The lower panels show the radii used for spectral 
extraction and the reduced x 2 of the spectral fits. 





Fig. 10. — Contours of the projected X-ray temperature are 
overlaid on the optical image of C1J1226. 9+3332. Contours are set 
at 6, 9, 12, and 15 keV. 



radius (kpc) 

Fig. 11. — The projected emission measure profile of 
C1J1226. 9+3332 as measured with Chandra and XMM-Newton is 
plotted along with the best fitting model. 



the high central temperature and lack of cool core in 
C1J1226. 9+3332, and is likely a result of the suspected 
merger activity in the system. 

The 3D gas density and temperature profiles enable the 
cooling time profile of C1J1226. 9+3332 to be computed, 
using the relation 



t 



cool 



1.5 x 10 10 yr 



.lO- 3 cm-V ViO 8 K ' (( " 
l)Sarazini ri986). The resulting profile is plotted in Fig 



T 



1151 along with the Hubble time at z = 0.89, and shows 
that the gas in the central regions of C1J1226. 9+3332 
would take ~ 6 Gyr to radiate its thermal energy away. 
This combined with possibility of a recent merger event 
(implied by the asymmetry in the temperature map and 
galaxy distribution) make the absence of cool gas in 
th e core of this system unsurprising. As pointed out 
bv lCagnoni et all $2001) there is a faint point-like radio 
source at the centre of C1J1226. 9+3332, consistent with 
a low-luminosity radio-loud AGN. This is also a candi- 
date for energy input into the cluster core at some point 
in the past. 
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Fig. 12. — Total (upper line) and gas mass (lower line) profiles 
of C1J1226. 9+3332. The data points indicate the midpoints of the 
temperature bins and the shaded regions and errorbars show the 
1<t confidence intervals from the Monte-Carlo randomisations of 
the data. The vertical lines indicates the radii of -R2500 an d i?soo- 



R (Mpc) 

Fig. 14.— The gas entropy of C1J1226.9+3332 is plotted as a 
function of radius with the shaded region indicating the la con- 
fidence intervals and the data points the measured values at the 
midpoint of the temperature bins. The dotted line shows the base- 
line entropy profile of C1J1226. 9+3332 due to gravitational effects 
alone IVoit c t al. 2005) and the solid line is the baseline profile 
with an additional constant entropy level (see text). 




R (Mpc) 

Fig. 13.— The gas mass fraction of C1J1226.9+3332 is plotted 
as a function of radius. The shaded region shows the la confidence 
interval and the vertical lines mark the radii of R2500 and Rsoo- 




R (Mpc) 

Fig. 15.— The cooling time of the gas in C1J1226. 9+3332 is 
plotted as a function of radius. The shaded region shows the la 
confidence interval and the horizontal line marks the age of the 
universe at 2 = 0.89. 



8. DISCUSSION 

The deep XMM-Newton observation of 
C1J1226. 9+3332 has enabled a uniquely detailed 
mass analysis of a z — 0.89 galaxy cluster. This allows 
comparison with theoretical density profiles for the first 
time at such high redshift. A profile was constructed 
using the densities calculated at the midpoint of each 
temperature profile bin, with uncertainties derived from 
the Monte-Carlo realisations. T hese data points were fit 
by a lNavarro etaD (|1996l 119971 NFW) profile: 



p(r) 



{r/r s )(l + r/r s y 



(7) 



where the normalisation S c can be expressed in terms of 
the halo concentration ca, and the overdensity factor of 



interest A: 



Sr = 



3 ln(l + c) -c/(l + c)' 



(8) 



The scale radius r s is then related to the overdensity ra- 
dius i?A by the concentration, ca = Ra/t s - This model 
has been used successfully to describe the density pro- 
files of both simulated and observed clusters of a range 
of masses. 

The model was fit to the density profile of 
C1J1226.9+3332 using A = 500, and the best-fitting pa- 
rameters were C500 
( x 2 / u = 4.2/4). Using A = 
tion to c 2 oo = 7.9+^4 



5.3±i;g and r s = 0.18 ± 0.04 Mpc 
200 scales the concentra- 
The data and best fitting model 
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Fig. 16.— The total and gas density profiles of C1J1226.9+3332 
are plotted, scaled by the critical density at z = 0.89 and -Rsoo- 
The solid line is the best fitting NFW profile to the data. 



Fig. 17.— C1J1226. 9+3332 is plotted on the V06 M - T relation. 
The masses have been scaled by the expected self-similar evolution. 



are plotted in Fig. EH The measured concentration of 
C1J1226. 9+3332 is significantly higher than that found 
for local clusters of a similar mass by V06. Furthermore, 
simulations predict that ca should decrease with increas- 
ing redshift (for a fi xe d mass) (e.g. iWechsler et al1l2f)f)2t 
iBullock et aLll200l . IWechsler et all <|2002|) also found 
that concentration increases with the time between the 
epoch at which a cluster is observed (z bs) and that at 
which it formed (z form ); ca oc (l+2 bs)/(l+2form)- With 
Zform ~ 2.6, the concentration of C1J1226. 9+3332 would 
be in line with these predictions. However, given the sus- 
pected merger activity in C1J1226. 9+3332, z b s ~ Zform 
seems more likely, with the high concentration resulting 
from the merger activity. 

8.1. The global properties of CU1226. 9+3332 and the 
scaling relations 

To enable comparisons with other works, we compute 
the global properties of C1J1226. 9+3332 within different 
radii. These are summarised in Table ^ The temper- 
atures and metal abundances were obtained from single 
temperature fits to the spectra extracted within each ra- 
dius using the combined XMM-Newton data. These rep- 
resent some of the most precise such measurements for 
a high-redshift system. The metal abundance is slightly 
lower tha n the values of Z k, 0.2,7 ,^ observed in local clus- 
ters (e.g. iFukazawa et aljfl998^ . but is consiste nt with 
the de clining abundance with redshift predicted bv lEttoril 
( 2005J) based on modelled supernovae rates. 

We computed the bolometric luminosity of the clus- 
ter by scaling the value mea sured within R^g n out to 
1.43 Mpc (for consistency with Markcvitch 1998) and out 
to infinity using the measured surface brightness profile. 
We find L x {r < 1.43 Mpc) = 4.57 ± 0.11 erg s" 1 and 
Lx{t < 00) = 5.12 ± 0.12 erg s _1 , including uncertain- 
ties on the temperature and metal abundance in those 
on the luminos ity. These luminosities are consistent with 
those found bv lMaughan et all l)2006|) based on a shorter 
XMM-Newton observation of C1J1226.9+3332. 

The measured properties of C1J1226. 9+3332 were then 
compared with local luminosity-temperature (L — T) re- 



lation of lMarkevitcbl (fl99 

L x (r < 1.43 Mpc) 
10 44 erg s" 1 



6.35E{z) 



kT \ 
i keV/ 



2.1)1 



(9) 



Here we have converted the L — T relation to our cos- 
mology, and included the predicted self-similar evolu- 
tion. The properties of C1J1226. 9+3332 were found to 
agree very well with the predictions of this relation. 
C1J1226. 9+3332 was then compared to the local mass- 
temperature (M — T) relation derived by V06 



1.61 



(10) 



The M — T relation is again converted to our cosmology 
with the predicted self similar evolution. The V06 data 
and M — T relation are plotted in Fig. El a l° n g with the 
data point for C1J1226. 9+3332, which falls significantly 
below the V06 M - T relation 

Finally. Ikravtsov et al.l l)2006|) recently proposed a new 
mass proxy for clusters, Yx, obtained by multiplying 
the gas mass and temperature measured within -R500 
together. This quantity was found to scale with mass 
with smaller scatter than other X-ray observables, and 
simulations showed that the scatter is less than 8% 
even for unrelaxed clusters, and that the scaling re- 
lation evolves self-similarly. The high quality X-ray 
data for C1J1226. 9+3332 enable us to test these pre- 
dictions observationally for the first time at high red- 
shift. In Fig. ITglwe re produce the Yx — M 50 o plot from 
iKravtsov et al.l 1120061) inc l uding the V06 data points. 
Following iKravtsov et all l)2006() we parameterise the 
Yx — M 50 o relation as 



E(z 



,2/5. 



M: 



10 14 M 



T.58 



o 



Fx 



4 x 10 13 M Q keV 



3/5 



(11) 



Here, the normalisati on found from the simulations of 
IKravtsov et al.l 1)2006(1 has been reduced by 15% (in line 
with the recommendation in that work), to better agree 
with the observed V06 cluster data. C1J1226. 9+3332 has 
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A r A (Mpc) kT(keV) Z(Z e ) M gas (10 la M Q ) A/ iot (lQi 4 M ) f g g S y x (lO 14 M0 keV) 
500 0.88 ±0.05 10.4 ±0.6 0.16 ±0.05 7-8to\' s 5 - 2 -o.s °- 15 ^o.oi 8A -o* 

2500 0.43 ±0.02 10.7 ±0.5 0.16 ±0.05 3.0 ± 0.2 3.2± u ;| 0.10 ±0.01 3.2 ±0.3 

TABLE 1 

Summary of the global properties of ClJ1226.9±3332 measured within different overdensity radii. The uncertainties 
include the statistical uncertainties and the effect of a 50% increase on the temperature errors to account for. 

calibration uncertainties. 
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Fig . 18.— C1J1226.9±3332 is plotted on the IKravtsov et al.l 
(2006) Yx — M500 relation along with the V06 data. The masses 
have been scaled by the expected self-similar evolution. 



a mass that is slightly low compared to the local clus- 
ters (although it is within the scatter at the la level). 
This could be due to stronger than expected evolution in 
the Yx — M500 relation, with E{z)" ji / b resulting in good 
agreement. A more likely explanation though, is that 
the offset from the relation is caused b y errors in one 
or bot h parameters. The simulations of IKravtsov et aD 
( 2006) show that the scatter in Yx — M 50 o is very small, 
even for merging clusters. Taking the measured Yx of 
C1J1226. 9+3332 as an indicator of the true total mass 
then implies that the X-ray analysis underestimates the 
true mass by ~ 30%. Such a difference would also im- 
prove the agreement of C1J1226.9+3332 with the M - T 
relation. 

Numerical simulations have shown that merging clus- 
ters tend to have lower temperatures for a given 
mass than relaxed systems, as the increase in ther- 
mal energy of the gas doe s not match the increase 
in mass due to the merger llMathiesen fc Evrardll200ll 
IKravtsov et alJl2"006j) . Thus the suspected merger event 
in C1J1226. 9+3332 could be responsible for an under- 
estimate of the total mass using X-ray methods. The 
good agreement of C1J1226.9+3332 with the L - T rela- 
tion is not inconsistent with this explanation, as merg- 
ers have been found to boost the luminosity and tem- 
perature of clusters approximately simultaneously, mov- 
ing them along the L — T relation (R andall et al J 120021: 
iRowlev et aIll2004D . 

8.2. The effect of systematic temperature errors 



As discussed in fl5.2l an artificially high absorbing col- 
umn was required when fitting the MOS spectra to pro- 
duce agreement with the PN data. It appears likely that 
the PN data are the more reliable in this case, as the 
best-fitting PN absorbing column is close to the Galac- 
tic value while the MOS is significantly higher. How- 
ever, the full analysis was also repeated with the MOS 
column fixed at the Galactic value, resulting in system- 
atically higher t emper atures in the combined EPIC fit. 
Vikhl inTn et alJ 1(20051 ) also found cluster temperatures 
measured with Chandra that were systematically higher 
by ~ 10% than those measured by XMM-Newton. It is 
thus instructive to evaluate the effect that systematic in- 
crease to the individual temperatures would have on the 
inferred properties of C1J1226. 9+3332. 

Reassuringly, the total mass, gas mass, overdensity 
radii, and luminosity of C1J1226. 9+3332 did not change 
significantly when the Galactic column was used for the 
MOS data. These quantities are all robust to residual un- 
certainties in the EPIC calibration. The only significant 
change was in the temperatures, with the global tem- 
perature increasing to 11.8 keV, with a corresponding 
increase in Yx- The agreement with the L — T, M — T, 
and Yx — M500 relations were thus all worsened by the 
higher kT. The derived mass did increase slightly in 
this test, but remained within our uncertainties because 
we included a large systematic contribution in our er- 
ror budget and because the shape of the temperature 
profile did not change significantly in this experiment, 
even though its normalisation did. An isothermal mass 
estimate would have been much more sensitive to the 
increase in kT. 

9. SUMMARY AND CONCLUSIONS 

High quality XMM-Newton and Chandra data of 
C1J1226. 9+3332 have enabled the most precise X-ray 
mass analysis of any such high-redshift system. The tem- 
perature profile of the system showed no sign of central 
cooling, with a hot core and a radially declining pro- 
file. The cooling time of the gas only drops marginally 
below the Hubble time in the centre of the cluster. A 
temperature map showed asymmetry in the temperature 
distribution with a hot region that is apparently associ- 
ated with a subclump of galaxies, but is not visible in 
the X-ray surface brightness. This is likely to be result 
of a merger event in the cluster, but does not appear to 
significantly affect the overall temperature profile. 

The entropy profile in C1J1226. 9+3332 was not well de- 
scribe d by the gravitational baseline profile of lVoit et alJ 
( 2005) and required an additional constant entropy level. 
The /g as profile increases with radii out to i?5oo, and 
is consistent with the profiles observed in local clusters. 
The global properties of C1J1226. 9+3332 were compared 
with the local scaling relations and found to be consis- 
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tent with the L — T relation, but to fall below the M — T 
and Yx — M500 relations. As Yjc is insensitive to merg- 
ers, this implies that the total mass is underestimated 
by ~ 30% in the X-ray analysis, probably due to de- 
partures of the gas from hydrostatic and virial equilib- 
rium. The high NFW concentration parameter found for 
C1J1226. 9+3332 is likely also due to the merger activ- 
ity. Given the high luminosity of C1J1226. 9+3332, it is 
likely to be the only distant cluster for which such rich 
X-ray data are obtained with the current generation of 
satellites. 
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APPENDIX 

A SIGNIFICANCE CONTOURING METHOD 

Here we describe the procedure used to define contour levels based on the significance of the emission enclosed between 
them. For this method to succeed, the image used to define the contour regions needs to be sufficiently smoothed that 
smooth contiguous regions can be defined even where the number of counts are low. An adaptive smoothing algorithm 
is most suited for this purpose as it uses a large smoothing kernel in regions of low signal to noise, but maintains high 
spatial resolution w here the signal to noise ratio is high. Two adaptive smoothing algorithms were tested; asmooth 
labeling et alJl200l upon which the CIAO task csmooth 6 is based), and the XMM-Newton Science Analysis Software 
(SAS) smoothing tool (this tool is also called asmooth, but for clarity we refer to it herein as xmmsmooth) . We 
found that for this particular purpose, the asmooth was more suitable than xmmsmooth as xmmsooth created artifacts 
around point sources and pixelation in the cluster cores (even with a high smoothing significance threshold) that were 
incompatible with our contouring algorithm. We note that both asmooth and xmmsooth are less sensitive to regions 
of depressed emission than regions of excess emission. 

To set the lowest contour level, the flux was measured in the background image for that cluster (after it had been 
normalised to match the local background level) in a circle of radius 6" at the cluster centre. The lowest contour level 
was then set at a flux of 3cr above this background level. The flux interval between this level (/&) and the peak of the 
cluster flux was divided into 500 logarithmically spaced intervals (/<; i = 1 — 500). The next contour level was set at 
fi , and spatial regions bounded by the contours at fa and f\ were defined on a smoothed image of the X-ray emission. 
For each distinct region, the significance of the emission contained in that region above the background was measured 
using the counts in the original unsmoothed source and background images. Regions that contained emission detected 
at > 3<7 significance were kept, and the others were rejected. 

The process was then repeated with an upper flux limit of /-j. Where regions had passed the previous step, new 
"child" regions were defined bounded by the contours at f\ to f 2 , and were by definition fully enclosed by the passing 
parent region from the previous step. Where regions failed the previous step, new regions bounded by contour levels 
fb and /2 were defined. This process was then repeated for all of the flux levels fi, with the significance of the emission 
enclosed by each region measured with respect to its local background. The local background for all child regions 
was measured in their respective parent region. For the lowest level regions bounded by contours at /{, and fi, the 
blank-sky background image was used to measure the background level. 

An i mage o f some of the regions used to define the contour levels in the Chandra image of C1J1226. 9+3332 is shown 
in Fig. IA19I for illustration. Region 1 is bounded by contour levels at fb and f\ , and contains emission detected at 
> 3(7 above the emission in region 1 in the background image. Region 2 uses the emission in region 1 as the background 
level above which to measure its significance. Regions 3 and 4 use region 2 as a background, but their significance 
is measured independently. Note that this contouring scheme means that contour levels cannot be interpreted in the 
same way as in standard contour plots; separate contours that are the same number of contours above the background 
are unlikely to correspond to the same flux level in the smoothed image. 

When applying this method to a large number of Chandra observations of clusters (Maughan et al. 2006 in prep.) it 
was found that an artificial "banding" could result. This is a repeated sequence of two tightly grouped contours followed 
by two widely spaced contours. The reason for this is that the statistical uncertainties on the blank-sky background 
counts in a region are generally very small, so a narrow contour band is sufficient for a 3<r detection. When that 
narrow band is used for a background for its child region, the uncertainties on the background level are relatively large 
due to the small number of counts, so a wide contour band is required for a 3c detection. The next band is then 
narrow because of the relatively small background uncertainties in the wide parent region, and the process repeats. 
In order to prevent this, an additional criterion was imposed to regulate the process. In order for a region to pass, 
it was required that the contained emission be detected at > 3<r and that the statistical uncertainty on the enclosed 
flux be smaller than some threshold. A threshold of 5% was used for the majority of the clusters, but a higher (lower) 
threshold was required for the faintest (brightest) clusters. This regulation prevents the banding problem and means 
that, particularly in the brightest clusters, the significance of some regions is somewhat larger than the 3cr threshold. 



6 N ote however that the CIAO csmooth is inferior to the Ebcline et al. (2006) asmooth algorithm. This is addressed by Ebcliim et al. 
(2006) and some shortcomings of csmooth are pointed out bv lDiehl fc Statlerf 12005). 
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